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We analyze the superconducting state and c-axis charge dynamics of cuprates using a charge or- 
dered bilayer superlattice model in which pairing is supported by inter-layer Coulomb energy gain 
(potential energy driven superconductivity). The superlattice nature of high-T c superconductivity 
is experimentally suggested by the smallness of the in-plane coherence length £ a f, ~ 10 — 30A which 
is comparable with a width of a 3ao x 3ao to 8ao x 8ao (ao ~ 3.9A) square supercell lattice layer. 
The 2D pair-condensate can be characterized by a charge ordered state with a "checkerboard" like 
pattern seen by scanning tunneling microscopy. The 2D 3D quantum phase transition of the 
hole-content at T c , supported by c-axis optical measurements, is also studied. The pair condensation 
might lead to the sharp decrease of the normal state c-axis anisotropy of the hole content and hence 
to the decrease of inter-layer dielectric screening. The drop of the c-axis dielectric screening can be 
the primary source of the condensation energy below T c . We find that a net gain in the electrostatic 
energy occurs along the c-axis, which is proportional to the measured condensation energy (Uo) and 
with T c : E^ D ~ 2N 2 Uo ~ ksT c and is due to inter-layer charge complementarity (charge asymmetry 
of the boson condensate) where N is the real space period of the condensate. The bilayer model 
naturally leads to the effective mass of m* m 4m e found by experiment. The static c-axis dielectric 
constant e c is calculated for various cuprates and compared with the available experimental data. 
We find correlation between T c and the inter-layer spacing d, e c and with the coherence area of the 
condensate. PACS numbers: 74.20.-z, 74.25.-q, 74.72.-h 



I. INTRODUCTION 

It is more or less generally accepted now that the con- 
ventional electron-phonon pairing mechanism cannot ex- 
plain cuprate superconductivity, because as high a tran- 
sition temperature as 16AK (the record T c up to now 
[1]) cannot be explained by the energy scale of lattice 
vibrations without leading to lattice instability [2]. It 
is already well established that much of the physics re- 
lated to high-temperature superconductivity (HTSC) is 
in 2D nature, one of the basic questions to be answered, 
however, in the future is whether HTSC is a strictly 2D 
phenomenon or should also be described by a 3D theory. 

A well known experimental fact is that the zero re- 
sistance occurs along the ab plane and the c-axis at the 
same critical temperature [3] suggesting that there must 
be inter-layer (IL) coupling involved in the mechanism 
which drives the system into HTSC. Another experi- 
ment, such as the intercalation on Bi2212 [4], however 
leads to the opposite conclusion. Intercalation of I or 
organic molecules, which expands the unit cell signifi- 
cantly along the c-axis, does not affect T c . This find- 
ing is against IL coupling and supports low-dimensional 
theories. The large anisotropy of the resistivity (and 
of other transport properties) is again not in favour of 
3D theories of HTSC [3]. Moreover, Basov et al. re- 
ported c-axis optical results and detected a small low- 
ering of kinetic energy in the IL transport of cuprates 
[5]. Experiments on Y Ba^CuzO^-y (YBCO) ultra- 



thin artificial HTSC compounds, sandwiched between 
thick nonsuperconducting PrBa^Cu-iO^-y (PBCO) lay- 
ers [3,6] and on (BaCuC>2+x)2/ ' {GaCuOi) n [7] indi- 
cate the continous decrease of T c with the decreas- 
ing number of superconducting (SC) layers. In par- 
ticular the one- unit-cell thick sample of YBCO/PBCO 
heterostructure exhibits no HTSC with the T c s» 
20 K [6,8]. Other heterostructures, such as the 
(Ba . 9 Ndo.iCu02+ x h/{CaCu02h/(Bao.9Ndo.iCu0 2+x ) 
containing a single bilayer SC block isolated from each 
other by insulating blocks, were shown to have T c w 55 
K [9] . These findings again indicate the importance of IL 
coupling in high temperature superconductors. There are 
a couple of other findings against and supporting the 3D 
nature of HTSC [10-12]. Most notably the systematic 
dependence of the transition temperature T c on the o 
axis structure and, in particular, on the number of CuO-2 
planes in multilayer blocks also strongly in favour of the 
3D character of HTSC. It is therefore, a fundamental 
question whether at least a weak IL coupling is needed 
for driving the system to a superconducting (SC) state 
or a single CuO^ layer is sufficient for HTSC [13]. 

There has been considerable effort spent on under- 
standing HTSC within the context of IL coupling mecha- 
nism in the last decades in which c-axis energy is available 
as a pairing mechanism [14-17]. In other approaches the 
importance of IL hopping is emphasized vs. the direct 
IL Coulomb interaction of charged sheets [10,13]. There 
is another theory providing explanation for HTSC using 
the general framework of BCS combined with IL coupling 



[18]. The so-called IL tunnelling (ILT) theory [10,13], 
however, is no longer considered a viable mechanism for 
SC in cuprates since ILT could provide no more then 1% 
of the condensation energy in certain cuprates [2,11,19]. 

Recently, optical data have been reported for Bi2212 
[20], which supports the quantitative predictions of 
Hirsch [21], in which the carriers lower their kinetic en- 
ergy upon pair formation (kinetic-energy-driven super- 
conductors) and hence this is the largest contribution 
to the condensation energy. Anyhow, the role of direct 
IL Coulomb interaction is still not ruled out as a pos- 
sible explanation for HTSC, especially if we take into 
account that no widely accepted theory is available until 
now which accounts for the apparent 3D nature of HTSC 
[10,11]. 

In this paper we propose a simple phenomenologi- 
cal model for explaining the 3D character of HTSC in 
cuprates supported by calculations. We would like to 
study the magnitude of direct Coulomb interaction be- 
tween charge ordered square superlattice layers as a pos- 
sible source of pairing interaction. Our intention is to 
understand HTSC within the context of an IL Coulomb- 
mediated mechanism. Of particular relevance to our in- 
vestigation are the doping, multilayer, and pressure de- 
pendence of T c in terms of a charge ordered superlattice 
nature of pair condensation, which is often neglected in 
the past. 

On the other hand the size of a characteristic super- 
lattice can directly be related to the in-plane coherence 
length £ a b of cuprates, which is proportional to the lin- 
ear size of the pair condensate (real space pairing) in the 
ab-plane [22] . The relevant length scale for superconduc- 
tors is the characteristic size of the Cooper-pair and can 
be estimated by un uncertainty-principle argument [22]. 
Only those charge carriers play a decisive role in HTSC 
which has energy within ~ ksT c of the Fermi energy and 
sets in at T c with the momentum range Ap w )cbT c /vf, 
where vf is the Fermi velocity, leading to the definition 
of the Pippard's characteristic length [22] 



where a is a numerical constant of order unity, to be 
determined. t; a b is a relevant number for describing the 
2D confinement of the pair-condensate wave function be- 
low T c . £ a b can directly be obtained according to the 
anisotropic Ginsburg-Landau theory via the measure- 
ment of the upper critical field H C 2, a b [22]. 

The IL charging energy we wish to calculate depends 
then on the IL spacing (d), the IL dielectric constant 
e c , the hole content p and the size of the superlattice. 
A large body of experimental data are collected which 
support our model. A consistent picture is emerged on 
the basis of the careful analysis of this data set. Finally 
we calculate the static c-axis dielectric constant e c for 
various cuprates which are compared with the available 
experimental observations. 



It is commonly accepted that charge carriers are 
mainly confined to the 2D CuOi layers and their con- 
centration is strongly influenced by the doping agent via 
hole doping. Holes (no charge and spin at a lattice site) 
in the 2D CuOi layers are the key superconducting el- 
ements in high temperature superconductivity (HTSC). 
A characteristic feature of many high temperature super- 
conductors (HTSCs) is the optimal hole content value of 
p w 0.16e per CuOi layer at optimal doping measured 
in the normal state (NS) [23]. This general feature of 
cuprates can be summarized in the parabolic dependence 
of maximum critical temperature T™, 

T c /T c m = l-82.6(p- Po ) 2 , (2) 

where TJ™ corresponds to the optimal hole concentration 
p [24]. Since T c appears to be maximized at p a w 0.16, 
we pay special attention to the IL Coulomb interaction 
in cuprates at this particular hole concentration. 

Our starting hypothesis is that the hole content goes 
through a reversible 2D 3D quantum phase transition 
at T c in layered cupper oxides. Optical studies on the c- 
axis charge dynamics reveals this phenomenon: the c-axis 
reflectance is nearly insulating in the NS but below T c is 
dominated by the Josephson-like plasma edge [5] . Below 
T c a sharp reflectivity edge is found at very low frequen- 
cies (lower then the superconducting gap) for a variety 
of cuprates arising from the carriers condensed in the 
SC state to the afe-plane and from the onset of coherent 
charge transport along the c-axis [25-29]. The appear- 
ance of plasma in the SC state and the absence of it in 
the NS seem to be a common feature of HTSC cuprates 
[28]. Band structure calculations predict an appreciable 
c-axis dispersion of bands close to the Fermi surface and 
thus an anisotropic three-dimensional metallic state [30] . 
Other first principles calculations indicate that above T c 
the hole charge p is charge transferred to the doping site 
(2D — ► 3D transition) [31,32]. Furthermore, we expect 
that below T c the hole-content condenses to the sheets 
forming anti-hole regions (hole-content charge at a lat- 
tice site, 3-D — ► 2-D transition). An anti-hole corre- 
sponds to an excess charge condensed to a hole lattice 
site in the sheet below T c . Therefore, in the normal state 
hole doping is the dominant charge transfer mechanism 
while in the superconducting (SC) state the opposite is 
true: the IL hole-charge is transferred back to the CuOi 
planes (charging of the sheets). The latter mechanism 
can be seen by the measurement of transport proper- 
ties as a function of temperature [3,33-35] if we assume 
that the increase in the density of in-plane free carriers 
below T c is due to pair condensation. The sharp temper- 
ature dependence of the c-axis dielectric constant e c and 
optical conductivity [36] seen in many cuprates and in 
other perovskite materials also raises the possibility of a 
2D 3D charge density condensation mechanism at T c 
[3,38]. Therefore, the pair condensation can be described 
by an anisotropic 3D condensation mechanism, and by a 
doping induced 2D-3D dimensional crossover [12]. 



In any naive model of electron pairing in cuprates the 
Coulomb repulsion is troublesome. When the pairs of 
charged carriers are confined to the 2D sheets, naturally 
a net self-repulsion of the pair condensate occurs. Al- 
though short-range Coulomb screening in the dielectric 
crystal can reduce the magnitude of the repulsion but it 
is insufficient to cancel Coulomb repulsion [39]. Spatial 
separation should reduce the interaction strength, but 
even at 14A e 2 /r leV if unscreened. One can as- 
sume capacitative effect between the CuOi planes: The 
charged boson condensate in one plane is stabilized by 
a deficiency of that charge on another plane [39]. The 
various forms of the capacitor model are associated with 
the 3D character of HTSC, considering the inter-layer 
charge reservoir as a dielectric medium [39]. The inter- 
layer charging energy might be insufficient to stabilize 
the self-repulsion of the holes in one plane and the self- 
repulsion of the charge condensate on another plane in 
the capacitor model. Also, the superconducting proper- 
ties of the hole-rich and nearly hole-free sheets would be 
different which is not verified by experimental techniques. 

Our intention is to combine the 2D and 3D nature of 
HTSC. Therefore, we propose a phenomenological model 
in which the charge distribution of the planes is polar- 
ized in such a way that holes and anti-holes (hole-electron 
pair) are phase separated within each of the sheets lead- 
ing to a charge ordered state (COS). Recently, signatures 
of charge ordering have been found in various cuprates 
and manganites in the presence and the absence of mag- 
netic fields [41-44]. Furthermore, the "checkerboard" 
charge pattern seen in Bi2212 [41] with spin periodicity 
of ~ 8ao strengthens our expectation that the SC state 
can be characterized by a COS. We focus on, then the 
magnitude of IL coupling (direct IL Coulomb interaction) 
between 2D static charge ordered superlattices. 

Charging of the sheets: The CuOi layers carry neg- 
ative charge obtained from the charge reservoir even in 
the insulating stochiometric materials, e.g. in the infi- 
nite layer compound CaCuOi each layer is charged by 
2e charge donated by the Ca atoms. In (CuOi) 2 ~ then 
the electron configuration of Cu and O is 3d 10 4s 1 and 
2p 2 z . Unit (Cu0 2 ) 2 ~ is typical of any undoped (stochio- 
metric) cuprates and ab initio calculations provide ap- 
proximately the charge state Cu 15+ Of 5 ~ [3,30] which 
is due to charge redistribution between Cu and O. The 
(CuOi) 2 ~ plane itself is antifcrromagnetic and insulat- 
ing, e.g. there are no holes in the "overcharged" layers. 
Upon e.g. oxygen doping, however, the doping charge p 
is transferred along the c-axis to the doping site, since 
the doping atom exhibits a rare gas electron structure 
(O 2 -) [31,32]. 

Our basic assumption is that in the SC state every sec- 
ond hole with the charge of +p is filled up by the doping 
charge p due to the condensation of the hole-content p 
to the sheets. Briefly the c-axis anisotropy of the hole 
content is strongly temperature dependent in cuprates. 
The phase separation of the holes and anti-holes (lead- 
ing to a charge ordered state) is stabilized by the intra- 
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FIG. 1. Upper panels: The charge odered state of the type 
of a "checkerboard" of the hole-anti- hole condensate. Opened 
and filled circles denote the holes with a charge of qi — +0.16e 
and anti-holes (qj — — 0.16e), respectively in the 4ao x 4ao 
(5 x 5) square lattice layer model. Note that the left panel 
accomodates 13 anti- holes which corresponds to 2e + 0.16e 
charge. The right panel contains 12 anti-holes correspond- 
ing to 2e — 0.16e charge. The two lower panels correspond 
to the antiferromagnetic insulating state (left) and to the 
holed-doped system (normal state, right). 




FIG. 2. The charge ordered state of the hole-anti-hole con- 
densate in the bilayer-hole 4ao x 4ao superlattice model. Note 
the charge asymmetry between the adjacent layers. The bi- 
layer can accomodate a pair of boson condensate (4e). Note- 
worthy that holes (empty circles) and anti-holes (filled circles) 
can be characterized by stripes along the diagonal lines. 



and inter-layer hole- anti-hole interactions. Basically the 
phase separation of holes and anti-holes in the planes is 
the manifestation of strong Coulomb correlation in the 
SC state. The superconducting charge-ordered states in 
cuprates are characterized theoretically based on the as- 
sumption that cuprates are in close proximity to quan- 
tum critical points where spin/or static charge order oc- 
curs in the SC state [45]. In this paper will focus on the 
Coulomb interaction between the charge ordered sheets 
screened by the IL dielectric media (charge reservoir) . It 
will be immediately apparent from our analysis that net 
energy gain occurs from the IL Coulomb interaction of 
charge ordered states when the hole-electron charge pat- 
tern is asymmetrically distributed in the adjacent layers 
(Fig. 2). ' 

II. THE SUPERLATTICE MODEL 

We propose to examine the following superlattice 
model of pair condensation: A pair of charge carriers 
(2e) is distributed over 2/0.16 = 12.5 CuC>2 unit cells 
in a square lattice layer if the 2e pair is composed of the 
hole charge p s=s O.I6/CW.O2. However, allowing the phase 
separation of hole-anti-hole pairs, every second unit cell 
is occupied by — 0.16e (anti-hole), and the rest is empty 
(holes, +0.16e), therefore we have 25 unit cells for a con- 
densed pair of charge carriers (Fig 1., that is the unit lat- 
tice of the pair condensate in the 5x5 supercell model) . 
Therefore, 2p D = — 0.32e hole charge condenses to ev- 
ery second hole forming anti-hole sites. The remarkable 
feature is that the size of the 5x5 condensate (four lat- 
tice spacings, 4ao ~ 15. 5A) is comparable with the mea- 
sured small coherence length £ b of single-layer cuprates 
(Cab ~ 10l to 20l) [3,22,44,47]. £ afc can directly be re- 
lated to the characteristic size of the wave-pocket of the 
local Cooper pair (coherence area) [3,22]. The supercell 
model can be applied not only for p a w 0.16e but also 
for the entire doping regime. Our expectation is that the 
COS of the 5x5 model given in Fig. 1 can be an effective 
model state for describing the SC state. An important 
feature of this model is the charge separation dq in the 
charge ordered state, where +0.16e and — 0.16e partial 
boson charges are localized alternatively (dq — 0.32e). 
The hopping of charge carriers from the anti-hole sites to 
the holes can reduce the magnitude of the charge sepa- 
ration dq leading to the extreme case when dq w which 
is nothing else then the (CuC>2) 2 ~ antiferromagnetic in- 
sulating state. Therefore characteristic quantity of the 
SC state dq is directly related to the hole-content seen in 
cuprates in the NS as a function of doping. 

The hopping of anti-holes in such a charge ordered lat- 
tice layer might lead to SC without considering any lat- 
tice vibration effects when dq is in the optimal regime. 
The freezing of the COS leads to insulating Wigner crys- 
tal (the NS) [40]. Therefore the melting of the crystal- 
lized COS is required for SC [45]. HTSC can be char- 



acterized in this way as a competition between liquid- 
crystal-like striped COS (SC) and Wigner-crystal-likc 
phases (NS) in a semiclassical theory of hole dynamics 
[40]. 

Below T c the charged-ordered state becomes stable 
compared with the competing phase of the NS supported 
by IL charging energy. Holes and anti-holes are placed 
in such a way in adjacent layers to maximize IL charging 
energy. Therefore, holes in one of the layers are always 
covered by anti-holes in the other layer and vice versa 
(FIG. 2, inter-layer electrostatic complementarity, bilayer 
5x5 (4ao x 4ao) model). An important feature is then 
that the boson condensate can be described by an inter- 
layer charge asymmetry. The IL coupling of boson-boson 
pairs in the bilayer 5x5 model naturally suggests the 
effective mass of charge carriers m* 4m e , as it was 
found by measurements [48,51]. In the next section we 
generalize the 5x5 model to represent a real space pe- 
riodicity of N x N coherence area. The charge ordered 
state presented in this article can also be studied as a 
charge density wave (CDW) condensed on a superlattice. 
Using the CDW terminology one has approximately two 
CDWs within a characteristic bilayer with an amplitude 
of ~ 0.16e. This system might not be insulating if we 
assume high rate of intersite hopping and low rate of 
IL hopping of the charge carriers. Then the charge or- 
dered state is stabilized due to the attractive IL Coulomb 
interactions of asymmetrically condensed neighbouring 
CDWs. 

We characterize the purely hole-doped NS as follows: 
In each sheet each CuOi units posses p = +0.16e charge. 
The hole charge p is charge transfered to the charge reser- 
voir (doping site) above T c . In this model one can naively 
expect that the self-repulsion of a hole rich sheet can be 
quite large. In order to reduce the hole-hole self-repulsion 
of the sheets, a hopping of <~ p charge between the ad- 
jacent hole sites would result in the reduction of the re- 
pulsion leading to a Wigner-crystal-like COS depicted on 
FIG 3. As is well known, the hopping of certain amount 
of charge (~ p) between the holes leads to metallic con- 
ductance which is typical of hole-doped cuprates. This 
is the normal state (NS) of the cuprates. Due to the 
strong c-axis anisotropy of the hole charge in the NS IL 
dielectric screening is large in this state because of the 
increase in the dynamic component of the IL dielectric 
permittivity (e c ) and therefore the IL interaction energy 
vanishes, E^ s rts 0. 

Our proposal is to understand HTSC in such a way 
that the pair-condensate is stabilized at higher tempera- 
ture by the IL-charging energy provided by the antisym- 
metry of the pair-condensate charge density between the 
adjacent layers. 

The c-axis dielectric constant of the SC state is reduced 
to the average value of the IL dielectric medium due to 
the pair condensation of the hole charge to the sheets. No 
IL screening of the hole charge occurs in this state which 
may well lead to energy gain . When the system goes 
to SC state, the conduction band loses its 3D anisotropy 



FIG. 3. The charge ordered state of the normal state. 
Opened and filled circles denote the holes with a charge of 
qi « +0.32e and anti-holes (qj w — 0.16e), respectively in the 
5x5 square lattice layer model. 



due to the 2D pair-condensation of the hole-content. 



III. THE TOTAL ENERGY OF THE PAIR 
CONDENSATE 

Our intention is to develop a simple working hypoth- 
esis in which the IL charge reservoir provides an av- 
erage dielectric background and the hole content fur- 
ther enhances IL dielectric screening when 2D — > 3D 
transition occurs (hole-doping, the NS). The dielectric 
plasma provided by the hole charge results in the dy- 
namic screening effect of Coulomb interaction which is 
typical of hole doped cuprates in the NS. In the oppo- 
site case (3D — ► 2D, pair condensation) the IL dielectric 
screening is nearly reduced to the average background 
value of the IL ion core spacer. This is what leads to 
IL energy gain. The underlying source of the condensa- 
tion energy is then the energy gain due to the lack of 
dynamic screening in the SC state. The possibility of di- 
rect IL hopping of Cooper pairs is not considered within 
this model as it was ruled out as important aspect of 
HTSC [11]. 

We start from a very general description of our model 
system using e.g. a Hamiltonian similar to that is given 
elsewhere [15] or which can also partly be seen in gen- 
eral text books [56]. We would like to describe then the 
2D 3D condensation of the hole-content using the 
Hamiltonian 
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where Hf D is the BCS-type Hamiltonian of the intra- 
layer condensate. 
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In Eq. (4) cj^ i is the creation operator for electrons 
in the zth layer, with linear momentum k within the 



layer and spin a, using the effective-mass approximation 
£fe = h 2 k 2 /2m*. V = — \V\ is the attractive inter-layer in- 
teraction (3D coupling) and is assumed to originate from 
some of the proposed mechanisms [10,14-16]. The at- 
tractive IL term Hf- 3 given by 
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The first term in Eq. (5) is the direct IL hopping, while 
the second and third terms describe IL coupling in a 
general form, t is very small in oxide superconductors 
[11,15]. Y denotes the IL coupling assisted by direct 
Coulomb interaction between charged layers (this is of 
particular interest in our model). W is the coupling 
constant and can arise through the Coulomb interac- 
tion causing inter-band transitions at the Fermi surface 
to some of the occupied or empty bands away from the 
Fermi level, with finite dispersion, however, along the c- 
axis. In this paper we consider only the second term in 
Eq. (5) as the source of the attractive interaction for 
pairing and neglect the rest of the inter-layer Hamilto- 
nian Hff (t«0,W«0). 

We calculate then the energy of the nearly 2D electron 
pair condensate in the CuOi plane for the bilayer super- 
cell problem. It is assumed that the pair condensate be- 
haves as a nearly free electron gas with afr-kinetic energy 
Et%„ and potential energy, which is mainly its in-plane 



J kin 

self-electrostatic energy (E^ b ) and the out-of-plane inter- 
layer interaction energy E^ L . For simplicity, the rest of 
the electron system is neglected completely. The exter- 
nal potential of the condensate is also excluded in this 
model, that is the lattice-condensate interaction, which 
is assumed to be negligible in the charged (Cu02) 2 ~ sys- 
tem (at least its contribution is negligible to the conden- 
sation energy). In other words the ionic background of 
the planes is screened by the core and valence electrons 
of the {Cu02) 2 ~ plane. The kinetic energy of the charge 
condensate is due to the hopping of the charge carriers 
between the CuOi sites within the sheets. We do not take 
into account the complications due to ionic heterogene- 
ity, nonpointlikc polarization, etc. The effective Coulomb 
interaction between spinless point-charges may be ap- 
proximated by the expression V e ff(r) = e 2 /(47reoe c r), 
where e c takes into account phcnomcnologically the di- 
electric screening effect of the IL dielectric medium and 
confined hole charge. This kind of a rough approxima- 
tion has widely used by several groups in the last decade 
[14,39,48,65]. The important feature is that the c-axis 
dielectric screening (e c ) is nearly reduced to the static 
value of the average background dielectric constant in 
the SC state. In accordance with this a sharp temper- 
ature dependence of e c (u>) is found in BSCCO by c-axis 
optical measurements at T c [36]. The static out-of-plane 
dielectric function can be obtained from the sum rule 



e c = ei(0) 



e 2 (w') 



(6) 



where £2(0;) is the dynamic component of e c (uj) [63]. 

Although, the r-dependence of the screened in-plane 
electrostatic interaction between the condensed charges 
qi and qj is not perfectly 1/r, we approximate it with 
the expression ~ e 2 /47reoe a 6?" as well, where screen- 
ing is taken into account implicitly via e a b- The lowest 
eigenvalue of H given in Eq. (3) is E = (%!\H\^), where 

r 2) — J2k 9{k)e lk ' ri e lk ' r2 being the wavefunction of 
the interacting pair and g(k) is a pair-correlation func- 
tion [22]. Then the energy of the lowest eigenstate of the 
condensate in the SC state is approximated by using a 
pure dielectric form for the potential energy 
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where 7i is the Planck constant, m* ~ 4m e is the effec- 
tive mass for the holes induced in the half-filled bands 
[48], n, m are layer indices, r^ and r|"' m ^ arc the intra- 
layer and inter-layer point charge distances, respectively. 

The most important components are the interactions 
(1 2) 

with r\j' (bilayer components), however one has to 
sum up for the IL interactions with terms rj"' m \ where 
m = [2,oo]. Note that only the interactions of various 
layers with the basel bilayer (FIG. 2) are considered along 
the c-axis in both directions (up and down). The in-plane 
electrostatic screening e a \, is completely separated from 
the out-of-plane dielectric screening (e c ). qt,qj are the 
partial point charges/atoms in the N x N superlattice 
model at optimal doping. 
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where factor 2 is due to the fact that every second CuOi 
site is occupied by anti-holes, and each of them consists 
of 3 atoms (qij w ±0.053e at maximal charge separa- 
tion). For the sake of simplicity it is assumed that the 
charges are equally distributed among Cu and O atoms 
within a C11O2 site. The number of the lattice sites in 
the characteristic superlattice is iV 2 and 



U « (N - l)ao 



(9) 



where ao ~ 3.88A is the ab lattice constant. Nh — 2e 
is the charge of the electron pair. The anti-hole charges 
qahoie mus t satisfy the charge sume rule within a charac- 
teristic bilayer over a coherence area ~ £ 2 b 



N 

E 



q ahole = 



(10) 



where q? hole represents the anti-hole point charges. Nat- 

urally the charge neutrality Y^ij (li + Qj ) = 
is also required. In this paper we study the lattice size 
N x N — 5, which we found nearly optimal for a vari- 
ety of cuprates. However, it can be interesting to study 
the variation of the " characteristic" lattice size in differ- 
ent cuprates as a function of various parameters (dop- 
ing, pressure etc.). The kinetic energy of the boson 
condensate arises from the hopping of anti-holes (hole- 
anti-hole exchange) between the adjacent sites within the 
sheets, against the electrostatic background of the rest of 
the hole-anti-hole system. Single holc-anti-hole exchange 
is forbidden since extraordinary repulsion occurs which 
leads to the redistribution of the entire hole-anti-hole 
charge pattern. The collective intersite anti-hole hop- 
ping results in the kinetic energy of the condensate. If 
we chose appropriately the in-plane dielectric constant 
[20,58]), we get the following result: 
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Consequently, no net in-plane energy gain is available 
for HTSC in this model. In other words we expect no 
role of pairing induced kinetic energy gain contrary to 
other theories [2,10,21] since the kinetic energy of the 
condensate is cancelled by the in-plane electrostatic en- 
ergy. This is basically the manifestation of the virial 
theorem (2T + V = 0) for the a6-plane and is reasonable 
to expect its validity when weak IL coupling is assumed, 
where Ef m = 2T and Ef = V. If the inter-sheet cou- 
pling is not that weak, then kinetic energy gain can occur 
in the planes (the violation of the in-plane virial theo- 
rem). It can therefore be the subject of further studies 
that both afe-plane kinetic energy and the screened c-axis 
IL Coulomb interaction can be the source of the conden- 
sation energy in the SC state. 



IV. THE CONDENSATION ENERGY IN THE 
CHARGE ORDERED BILAYER-HOLE MODEL 

Important consequence of the model outlined in this 
paper is that the condensation energy (E con d) of the SC 
state can be calculated. Within our model the primary 
source of E con d is the net energy gain in IL Coulomb 
energy occurs due to the asymmetrical distribution of the 
condensed hole-charge in the adjacent layers (see Fig. 2) 
below T c . 

Per definition E con d is the free energy difference of the 
normal state and the superconducting state [11]. There- 
fore, taking the energy difference AE to t — \E^ a f — Ef£\ 1 
using Eq. (7) and dropping the terms coming from the 
a6-plane according to Eq. (11), the IL energy gain in the 
SC state (~ condensation energy) can be given as follows, 
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IL,SC\ 



(12) 



We can further simplify Eq. (12) if the NS contribution 
to Eq. (12) is El L ' NS « 0, which holds if IL coupling is 
screened effectively in the NS (large density of the hole 
content in the IL space, large e c (NS)). Then we have 



2N 2 U Q « E c 



E, 



IL.SC 



(13) 



where E\ L > sc is the Coulomb energy gain in the SC state. 
Uo is the experimental condensation energy given per 
unit cell. E con d is the condensation energy of the bilayer 
system with 2N 2 unit cells. 

The energy gain in the SC state is provided by the 
change in the inter-layer charging energy which is essen- 
tially the change of the out-of-plane hole-anti-hole inter- 
action energy below T c . The hole-conductivity in the SC 
state is strictly 2D phenomenon, no direct IL hopping 
of quasiparticles is considered within this model. T c is 
mainly determined by the inter-plane distance and by 
the static c-axis dielectric constant (the c-axis compo- 
nent of the dielectric tensor). An interesting feature of 
the N x N bilayer model is then that it is capable of re- 
taining the 2D character of superconductivity while T c is 
enhanced by 3D Coulomb interactions. 



V. THE CONDENSATION ENERGY AND T c 

There are number of evidences are available which 
suggest that HTSC occurs beyond the BCS limit. In 
those materials which contain nearly isolated single lay- 
ers, such as Bi2Ba 2 CuOe + s (T c = 20 K [64]) or su- 
perlattice structures (periodic artificially layered materi- 
als) such as O-doped {BaCuOi)2l(CaCu02) n thin films, 
when n = 1 [7] and in Y BaiCu-sO-i-s j PrBaiCuzOi-f, 
(YBCO/PBCO) [3,6,8] the critical temperature is lim- 
ited to T c < iQK (below the BCS limit), therefore 
these materials are not considered as high-T c super- 
conductors in this article. In these compounds the IL 
distance is so large that the CuC>2 planes are decou- 
pled and the superconducting properties can be under- 
stood within the BCS formalism. However, the multi- 
layer Bi-compounds (with the same charge reservoir), 
the {BaCu0 2 )2/{CaCu0 2 )n thin films, when n > 2 [7] 
or YBCO with thin PrBa 2 Cu 3 7 ^ s layer [3,6,8] exhibit 
HTSC. Therefore it is worth to explain the enhancement 
of T c beyond the BCS limit assuming other mechanism 
than the electron-phonon coupling. Inter-layer Coulomb 
coupling can be a natural source of the condensation en- 
ergy and of HTSC. 

Checkerboard-like charge pattern (COS) seen experi- 
mentally [41-43] directly leads to IL energy gain and to 
potential energy driven condensation energy if the hole- 
anti hole charge pattern is asymmetrycally condensed to 
the adjacent layers (Fig 2). Assuming that thermical 
equilibrium occurs at T c for the competing charge or- 
dered phases of the NS and the SC state the following 
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FIG. 4. The critical temperature (K) as a function of the 
bilayer condensation energy (2N 2 Uo, meV) using Eq. (14). 
The real space period N is directly related to £ a b via Eq. (9). 
The straight line is a linear fit to the data. 



equation for the condensation energy can be formulated 
using Eq. (13), 



2N 2 U Q 



Cab 
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Uo w k B T c a E, 



IL.SC 



(14) 



At a first look this formula seems to be unusual be- 
cause of the dependence of the condensation energy on 
T c . The available measurements of the condensation en- 
ergy on various cuprates show no correlation of Uo with 
the critical temperature. One of the important goals of 
this paper, however, to show that correlation can indeed 
be found with T c if 2N 2 Uq (the bilayer condensation en- 
ergy) compared with. In other words the condensation 
energy of the coherence bilayer-hole system shows corre- 
lation with T c . In order to test the validity of Eq. (14) 
we estimate the real-space period of the pair condensate 
using Eqs. (14) and (9), 
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Cab 
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(15) 



The results are given in Table I. N can directly be com- 
pared with the measured coherence length Cab via Eq. (9). 

Eq. (14) is clearly clarified in FIG 4 using mostly ex- 
perimental information for Cab, T c and Uo (values are 
given in Table L). For T12201 no measured Cab is found 
in the literature, therefore the estimated real space pe- 
riod is used (given in Table I.). The slope of the linear 
fit in FIG 4 precisely gives us ks which nicely confirm 
Eq. (14). The underlying physics of HTSC seems to be 
reflected by Eq. (14): the equation couples the variables 
T c , Uo and Cab- 

We find correlation between the condensation energy 
Uo and N: the larger Uo is connected with smaller 
N (see also FIG 5). The stronger localization of the 
pair- condensate wave function seems to lead to larger 
IL condensation energy and hence to larger T c . This 
is again an unexpected result, since the stronger local- 
ization of the Cooper pairs should increase the Coulomb 
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FIG. 5. The measured condensation energy (given in Table 
L, jueV/u.c.) as a function of the calculated coherence length 
(oo) using Eq. (15). 



self-repulsion of the condensate and hence should sup- 
press T c . Note, however, that within a liquid-crystal- like 
COS the self-repulsion problem is not crucial. The hole- 
anti-hole Coulomb interactions are attractive, although 
cancelled by the intra-hole and intra-anti-hole repulsions 
(intra CuO-i site repulsion, ) and by the in-plane kinetic 
energy via Eq. (11). 

The Uq value for LSCO is measured for a slightly un- 
derdoped sample [53] therefore a somewhat higher Uq 
might be obtained for the optimally doped material and 
hence smaller N and f ao can be calculated. Interestingly 
the most localized coherence area (~ 4ao) is provided 
by YBCO which is the less anisotropic material among 
HTSC cuprates (the resisitivity ratio p a b/p c ~ 200 [3]). 
The more anistropic Hgl201 gives us somewhat weaker 
localization of the condensate wave-pocket for similar T c . 
The comparison between various materials , however, is 
much more complicated. In general we can say that T c is 
a function of the following parameters considered in this 
study: N, p, the IL distance d and e c . In section VIII. we 
further analyse this complex behaviour of T c focusing on 
the calculated IL dielectric constant e c . The calculated 
real space periods given in Table I of the bilayer COS are 
in close agreement with the experimental coherence areas 
fab which supports the validity of our basic Eq. (14). 

The expression Eq. (14) leads to the very simple for- 
mula for the critical temperature 
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where Ni is the number of layers along the c-axis. When 
Ni — > oo, bulk T c is calculated. T c can also be calculated 
for thin films when TV; is finite and e c can also be derived 
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TABLE I. The calculated real-space period of the pair con- 
densate using the experimental condensation energies of var- 
ious cuprates and Eq. (15). 





T c (K) k B T c (meV) Uo (fieV/u.c.) 


N 


£ab(ao) 


LSCO 


32 


2.5 


13 a 


~ 10 


5-8 e 


T12201 


85 


7 


80 a 


~ 7 




Hgl201 


95 


7.8 


103 - 136 6 > 


-5-6 


5 b 


YBCO 


92 


7.5 


246 c 




3-4 / 


Bi2212 


89 


7.3 


107 d 


~ 6 


4-8 9 



where a c-axis average of e c is computed when Ni — > oo. 



Uo is the measured condensation energy of various cuprates. 
a from [53], b from [54] and £ ab from [55], c from [49], d from 
[24], e from [77], / from [22,50], 9 from [41,75], N is calculated 
according to Eq. (15) and £ a b is the in-plane coherence length 
given in ao w 3.88A. The notations are as follows for the com- 
pounds: LSCO (Lai.ssSro.iBCuO*), T12201 {Tl 2 Ba 2 Cu0 6 ), 
Hgl201 (HgBa 2 CuO i+ &), YBCO (YBa 2 Cu0 7 ) and Bi2212 
is Bi 2 Sr 2 CaCu 2 Os+s- 



VI. RELATIONS TO THE SUPERCONDUCTING 
GAP AND PAIRING 

Not useless to note again the correlation in Table I. 
between f ao and Uq. There seems to be a correlation 
between the real-space localization of the Cooper wave- 
function and the SC energy gain Uo. The real-space 
"shrinking" of the Cooper wave- function must directly 
be related to the enhancement of pairing. Within our 
model the primary source of increased pairing is the IL 
energy gain (Coulomb induced pairing). 

Optimally doped YBCO/PBCO with decoupled lay- 
ers show no HTSC (T c « 20 K) indicating that optimal 
carrier concentration within isolated layers leads to BCS 
superconductivity with low critical temperature. Setting 
in IL coupling (decreasing the thickness of the PBCO 
phase) T c is enhanced due to the increase of IL coupling 
[3,6]. Further experimental studies on artificially layered 
materials, such as the measuerement of f a b as a function 
of the thickness of the insulating phase PBCO should 
explain the importance of IL coupling in cuprates. The 
considerable increase in f ao as a function of IL decou- 
pling (increase in PBCO thickness or decrease in YBCO 
thickness) would be a strong evidence for the model pre- 
sented in sections IV. -VI. Searching for such an article in 
the literature an interesting publication is found which 
reveals our expectation [68]. Extensive measurements of 
the upper critical field in {YBCO/PBCO) n superlattices 
indeed results in the increase of ^ a b as a function of de- 
creased thickness of YBCO layers. f aD varies from 4 — 5ao 
up to ~ 8ao when the thickness of the YBCO phase is 
decreased from 10 unit cells to 3. Unfortunatelly this ar- 
ticle contains no measurements for 1 — 2 unit cell thick 
samples. 

In FIG 6 we give the measured coherence length in 
YBCO/PBCO getting the data from ref. [68]. We predict 
the rapid increase of f ao for dielectrically isolated single 
layers of YBCO for one or two-unit-cell thick samples 
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FIG. 6. The coherence length f a (, (A) as a function of the 
YBCO thickness (in unit cells) at fixed thickness of the barrier 
(insulating) layer PBCO (~ 48 and 96A) reproduced from 
Table II. of ref. [681. 



which is the reminiscent of purely BCS features. These 
striking results represent a strong evidence for the corre- 
lation between IL coupling and the coherence length £ a fc 
which is predicted by the model presented here. Strink- 
ingly, the contraction of the Cooper wave-function (the 
strengthening of pairing) is strongly coupled to the en- 
hancement of IL Coulomb coupling. 

As far as the microscopic nature of pairing is con- 
cerned, its relation to the SC gap should also be dis- 
cussed. In the BCS theory the SC gap is related to 
the pair binding energy of the Cooper pairs [22]. How- 
ever, keeping this scenario one can go easily into con- 
tradiction in HTSC. The increased gap ratio of cuprates 
2A w 8k B T c [3,67] can not easily be understood as the 
enhancement of pair binding energy (" pairing glue" ) . 
Within our understanding of HTSC the more or less 
generic gap ratio 2A ~ 7 — 8k B T c [3,67] is rather a char- 
acteristic feature of the superconducting COS and does 
not directly correlate with the pair binding energy in the 
Cooper-pairs. The recently found doping-induced linear 
shift of the gap ratio in Bi2212 [24,67,69] supports the 
lack of extra "pairing glue" in HTSC cuprates. The gap 
ratio monotonically decreases with doping showing not 
any sign of critical properties at the optimal or critical 
doping level. If the enhancement of pairing played a de- 
cisive role in HTSC then the gap ratio should peak at T™ 
which is obviously not the case [70] . The more or less con- 
stant value of the maximum gap ratio 2A ~ 7 — 9k B T c in 
various cuprates [3] is also not consistent with the " pair- 
ing glue" scenario. 2A should vary system by system as a 
function of T c if pairing enhanced with the increase of T c . 
However, systems with low T c w 40 K, such as the proto- 
typical HTSC material La2-x Sr x CuOi (LSCO) provides 
a relatively large gap as well (2A ~ 7.7k B T c ) [3,69]. Oth- 
ers measure a relatively small gap of 2A ~ 5k B T c ) for 
YBCO [72]. Anyhow, the interpretation of the measured 
gaps and its relation to pairing is still a matter of debate 
in HTSC [71,73]. We also attribute the pseudogap phe- 
nomenon of cuprates to the characteristic feature of the 



Wigner crystal-like COS of the normal state. This field 
is, however, obviously beyond the scope of the present 
study. 



VII. THE DIELECTRIC CONSTANT e c WITHIN 
THE BCS LIMIT 

The model outlined in the previous section allows us 
to estimate the static dielectric constant e c for various 
cuprates. First we give the brief derivation of the formula 
starting from the BCS equation for the gap [22], 



k B csk B T c = A(0), 



(18) 



where k B is the Boltzman constant and the gap ra- 
tio k B cs ~ 3.53 for normal superconductors at weak- 
coupling limit [22]. The BCS estimate of the conden- 
sation energy (the energy gain in the SC state) is pro- 
portional, however to ~ A 2 , therefore we can use the 
following expression [22], 



Econd = -^iV(Q)A(0) 2 = -\N(S))k 2 BCS k 2 B Tl 



(19) 



where iV(0) is the density of states at the Fermi surface 
[22]. We recall now Eq. (13), 
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Using Eq. (19) E con d can be substituted into Eq. (21) 
leading to a new expression, 
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N(0) can be obtained from first-principles calculations 
or from specific-heat measurements. N(0) is connected 
with the Sommerfeld constant 7 in the low-temperature 
electronic specific heat by the relation 



N(0) =7[2y*|(l + A)]" 1 , 



(23) 



A is a coupling constant due to strong electron-phonon 
renormalization. 7 can be obtained from specific heat 
measurements [3] 



1.437 = AC/T C 



(24) 



where 7 is usually given in [mJ/mol.K 2 ] therefore ~ 
2iV 2 7 must be used for the coherence area of ~ £, 2 b (twice 



for the bilayer). Finally we get the expression for the 
static c-axis dielectric constant 
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The only uknown parameter in Eq. (25) is A which is 
in the range of A = 1 — 2 when strong electron-phonon 
coupling is assumed [3]. 

We use later on in the discussion section this expression 
to calculate e c within the BCS limit for various cuprates 
comparing with the experimental dielectric constants. 



VIII. THE CALCULATED DIELECTRIC 
CONSTANT 

The calculation of the c-axis dielectric constants e c 
might provide further evidences for Eq. (14) when com- 
pared with the measured values [36,58,65]. In Table II 
we have calculated the static dielectric function e c using 
Eqs. (25) and (17) and compared with the experimental 
impedance measurements [59]. e c can also be extracted 
from the c-axis optical measurements using the relation 
[19,49] 
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where c, lo p and A c are the speed of light, c-axis plasma 
frequency and the c-axis penetration depth. 

In other cases only the e c of certain elements of the ion- 
core spacer is available such as e.g. BaO in Hgl201 and 
in YBCO. For the prototypical cuprate LSCO we get the 
value of e c = 27.3 which is comparable with the experi- 
mental value of 23 [58] using N — 7 which corresponds 
to 6a w 23.41 (£ ab w 20 - 30 A [77]). 

The overall good agreement of the calculated e c with 
the measurements is due to our finding that the IL charg- 
ing energy is surprisingly in the range of ksT c when the 
real space period N is reasonably chosen. Although we 
have no clear cut evidence for Eq. (14) we do not think 
that the agreement of various calculated properties with 
the measurements is accidental. In the rest of the article 
we collect further evidences which validates Eq. (14). 

Multilayer and pressure effects on T c can also be dis- 
cussed in terms of in-plane and out-of-plane localization 
of the hole charge. In multilayer systems E* L is com- 
posed of intra- and inter-block contributions. We use 
the notation block for the multilayer parts CaCuiO±, 
Ca-iCuzO^, etc. Eq. (13) can then be modified for mul- 
tilayer cuprates as follows, 
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TABLE II. The 
Eqs. (25) and (17) 
real-space period N 



calculated dielectric constant e c using 
in various cuprates as a function of the 
of the charge ordered state. 
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where N is the real space period of the characteristic square 
lattice (N ~ £ a f>/cio + l), the bold faced values are those which 
account the best for comparison with experiment and are in 
accordance with the results of Table I. dq is the charge 
separation in the charge ordered state, d is the inter-layer 
distance in A [14], T c is the experimental critical tempera- 
ture. e c (A) is from Eq. (25), e c (B) is from Eq. (17). b 
Cao.2,Sro.7Cu02 [61], c [58], or from reflectivity measure- 
ments using Eq. (26), uj p w 55cm _1 [79], A c « 3/wn [62], d 
experimental e e * v values are taken from Am Inst, of Phys. 
Handbook, McGraw-Hill, Ed. D, E. Gray (1982), the e c value 
of the ionic-background (Hgl201: BaO, T12201: TI2O5), The 
pressure dependent T c values and IL distance date are taken 
from [66,74], e from [19], 3 from [36], h for YBCO the reduced 
IL d = 8.5A is used instead of the c-axis lattice constant 
(inter-bilayer block distance), * from [65], 1 from [80], 3 from 
reflectivity measurements: lu v ~ 60cm~ 1 [79], A c w 0.9/im 
[37], e c can be deduced using Eq. (26). The notations are 
as follows for the compounds: LSCO (Lai.ssSro.isCuO^, 
Hgl201 (HgBa 2 CuO i+s ), T12201 {Tl 2 Ba 2 Cu0 6 ), Bi2201 
(Bi 2 Sr 2 Cu0 6 ) and YBCO (YBa 2 Cu0 7 ). 
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FIG. 7. The critical temperature T c (K, Eqs. (14)-(17)) vs. 
the number of unit cells along the c-axis in YBCO using the 
4x4 model. Circles and squares correspond to the calculated 
and experimental values [6]. 



1 denotes the number of layers. 

Therefore, intra-block charging energy further en- 
hances T c on top of the corresponding single-layer 
inter-block value. YBCO is a peculiar example of 
cuprates in which HTSC is purely coming from inter- 
block coupling as it was explained on the basis of 
Y Bci2Cu30r-s / PrBa,2Cu3C>7-s superlattice structures 
[3]. In YBCO the Y(Cu0 2 )2 2 bilayer alone does not 
show HTSC when isolated from each other in artificially 
layered materials (T c 2QK). The dependence of T c in 
YBCO on the number of unit cells along the c-axis is 
calculated using Eq. (16) and the results are depicted on 
FIG 7. The experimental points [6] are relatively well 
reproduced indicating that ~ 5 unit cell thick thin film is 
already the reminiscent of the bulk properties of cuprates. 

It is worth mentioning the intercalation experiment on 
Bi2212 [4]. Intercalation of I or organic molecules into the 
(BiO)2(SrO)2 layers of Bi2212 results in no significant 
change in T c . Basically inter-layer intercalation is intro- 
duced to reduce interlayer coupling (test of inter- layer 
theory, ILT) and to enhance anisotropy in SC properties. 
Using film deposition technique, ultrathin films of Bi-Sr- 
Ca-Cu-0 (Bi2212) have been synthesized [78]. The few- 
unit-cell-thick samples show HTSC similar to that of the 
bulk material independently of the film thickness. These 
results present great challenge to ILT [11] and support 
low-dimensional SC theories. 

In the light of the N x N model, however, it is not 
surprising that in Bi2212 the nearly isolated bilayers are 
are superconductors. The single layer material Bi2201 
{Bi 2 Sr2Cu0 6 ) gives very low T c (w 20K) [64], which 
indicates that the dielectrics {BiO)2{SrO)2 strongly re- 
duces IL coupling, indeed the estimated e c is quite large 
(Table II) which is attributed to the weakly interact- 
ing (BiO)2 bilayer (the BiO-BiO distance is ~ 3.7A 
[4]). Dielectric constant measurements [57] and the ob- 
tained small plasma frequency using c-axis optical mea- 
surements (u> ~ 5cm -1 ) [79] in Bi-compounds also pro- 
vide relatively large dielectric constants in accordance 



with our calculations. The extremely large conduction 
anisotropy 7 ~ 10 6 found in Bi-based cuprates also 
supports these findings [3]. Nearest-neighbouring Cu-0 
planes in Bi-compounds are nearly insulated [65] . There- 
fore, in the multilayer Bi-compounds bilayer and trilayer 
blocks are responsible for the high-T c , inter-block cou- 
pling is negligible. The multilayer Bi-compounds pro- 
vide then an example for pure intra-block HTSC. In these 
materials the multilayer-blocks are dielectrically isolated. 
However, in the most of the cuprates T c is enhanced both 
via intra- and inter-block effects. An important conse- 
quence of the bilayer N x N model is that isolated sin- 
gle layers do not show HTSC, coupling to next-nearest 
Cu0 2 plane is essential in HTSC. It is possible then to es- 
timate e c for the bilayer block, since E I c L ' mter ss in the 
Bi-compounds in Eq. (27). In the hypothetical bilayer 
system in Table II (Ca(Cu0 2 ) 2 , d = 3.191, T c = 89K), 
which is the building block of bilayer cuprates, we find 
the value of TV = 8 accounting for realistic e c . N = 8 
is in accordance with the measured relatively "large" co- 
herence length of £ a6 w 272 ((N - l)a = 7a » 27.3A) 
[75]. 

Also, upon pressure (p) the inter-layer spacing de- 
creases, which increases E^ L without the increase of the 
hole content. Saturation of T c (p) is reached simply when 
increasing IL sterical repulsion starts to destabilize the 
system. In systems, such as YBCO or LSCO, negative 
or no pressure dependece of T c is found [3] due to the 
short Cu-0 apical distance (dcuO ~ 2.4A) which leads 
to already steric repulsion at ambient pressure and to 
the weakening of HTSC. In these systems the net gain in 
charging energy is not enough to overcome steric repul- 
sions at high pressures. 

With these results, we are now in a position to reach 
the conclusion that the N x N model can readily ac- 
count for at least certain physical properties of multilayer 
cuprates, such as pressure, doping and multilayer depen- 
dence of T c . Furthermore it is possible to make some es- 
timations on the upper limit of T c using Eq. (17) for T c . 
Assuming relatively small e c ~ 10 and short inter-layer 
spacing d w 7.0A we get the value of T c w 333K for a 
strongly localized electron pair with 5x5 coherence area. 
Of course, we have no clear cut knowledge at this moment 
on how the pair condensate wave-function spreads upon 
varying e c and d. Effective localization of the Cooper- 
pair wave-pocket could lead, however, to room tempera- 
ture superconductivity under proper structural and dielec- 
tric conditions if the model presented above is applicable. 

According to the stripe scenario, the charge ordered 
state of the coherence area can also be seen in many 
cuprates as one-dimensional stripe order or charge den- 
sity waves [43]. The striped antiferromagnetic order 
found in LSCO by magnetic neutron scattering experi- 
ments [42] also implies a spin-density of periodicity ~ 8a 
a reminiscent of the coherence area. The incommensurate 
"checkerboard" patterns seen with a spatial periodicity 
of ~ 8ao in the vortex core of Bi2212 obtained by scan- 
ning tunneling microscopy [41] is also consistent with our 



N x N hole-anti-hole charge ordered state where N = 8 
to 9 in BSSCO. 

Finally we mention the recent results of Bozovic et al. 
[80] obtained for LSCO thin films under epitaxial strain. 
They reached the record T c = 51. 5K for 15-unit-cell thick 
film of LSCO on LaSrAlO^ substrate. The small vari- 
ation of the ab- or c-axis lattice constants in our model 
accounts for only 1 — 2 K increase in T c . We explain the 
more then 10 K enhancement of T c with the decrease of e c 
(Table II) under the conditions they used (O3 annealing, 
epitaxial strain provided by the substrate). 

IX. CONCLUSION 

In this paper we studied the pair condensation and 
confinement of the hole-content on a C11O2 superlattice 
layer as a function of inter-layer distance and dielectric 
permittivity of the charge reservoir. 

• The assumption of a 2D <=> 3D quantum phase 
transition of the hole-content at T c in HTSC mate- 
rials is thought to be an important general feature 
of pair-condensation and is supported by c-axis op- 
tical measurements and by first-principles calcula- 
tions. Our proposal is that the c-axis charge dy- 
namics of the hole-content contributes significantly 
to the condensation energy below T c . We find that 
the inter-layer capacitance is temperature depen- 
dent in cupratcs and therefore the drop of the c- 
axis dielectric constant can be seen below T c . This 
is what leads to then the stabilization of the super- 
conducting state vs. the normal state. 

• We have found that a pair condensate can be dis- 
tributed on a N x N square lattice layer in such 
a way that the lattice sites are filled by ±q = 
[0,0. 16e] condensed charge alternatively depend- 
ing on the hole(+g)-anti-hole(— q) charge separa- 
tion dq. In this way a charge ordered state of 
the pair-condensate occurs with a "checkerboard" 
like pattern seen recently by experiment [41]. The 
phase separation of hole-electron pairs (hole-anti- 
hole pairs) in this model is stabilized electrostat- 
ically. The maximum charge separation is dq w 
0.32e, if the optimal hole content p Q m 0.16e. 

• In the adjacent layer the electron-hole pairs are dis- 
tributed in a complementer way (charge asymme- 
try) in order to maximize the inter-layer charging 
energy. Holes on one plane are covered by anti- 
holes on the other, and vice versa. In this way 
we derived a charge ordered N x N bilayer-hole 
model of the superconducting state with inter-layer 
charge antisymmetry which directly leads to inter- 
layer Coulomb energy gain in the superconducting 
state. The IL charge asymmetry could directly 
be tested experimentally using e.g. two-cell-thick 



thin SC films, and the "checkerboard" STM images 
of both sides of the thin film could be measured. 
The complementer "checkerboard" charge pattern 
would reveal our expectations. 

• The superlattice nature of the pair condensate is di- 
rectly related to the smallest size of the condensate 
wave-pocket which is remarkably comparable with 
the measured in-plane coherence length of £ a & = 
10 — 20A (4a « 15. 6A) of single cuprates, where 
a w 3.9A is the in-plane lattice constant. The 
coherence area of the Cooper-pair wave- function in 
cuprates is strongly localized, which is due to inter- 
layer charging effects. 

• The bilayer model with 4e boson charge naturally 
implies the mass enhancement of m* w 4m e in ac- 
cordance with measurements. The calculated inter- 
layer charging energy is in the range of the ex- 
perimental condensation energy for the bilayer-hole 
system. 

• The static dielectric constant e c is calculated for a 
couple of cuprates and compared with the available 
experimental measurements. The general agree- 
ment is quite good indicating that the pure inter- 
layer electrostatic model leads to proper description 
of the static dielectric response of these layered ma- 
terials. 

• The basic microscopic mechanism of HTSC is to 
be understood within the BCS-Eliashberg the- 
ory. The limiting critical temperature for BCS- 
type superconductor is around 2QK as it was 
found for cuprates with nearly isolated CuOi layers 
(BSCCO, YBCO/PBCO superlattices etc.). The 
detailed study of this model showed that the inter- 
layer charging energy is proportional to the ther- 
mal motion at T c , if the c-axis dielectric constant 
e c and the coherence area is appropriately chosen. 
The correlation between the coherence area £ a b and 
inter-layer coupling is predicted by our model. The 
stronger IL coupling leads to smaller This re- 
lation is evidenced by the measurements of £ a {, in 
YBCO/PBCO films with varying YBCO thickness. 

If the physical picture derived from our model is cor- 
rect, it should be a guide for further experimental studies 
aiming to improve SC in cupratcs or in other materi- 
als. This can be done by tuning the IL distance and e c 
(increasing the polarizability of the dielectric, hence de- 
creasing e c ) in these materials. The application of this 
model to other class of HTSC materials, such as fullerides 
or MgB 2 is expected to be also effective. 
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